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ABSTRACT 

Supernova explosions inject a considerable amount of energy into the interstellar 
medium (ISM) in regions with high to moderate star formation rates. In order to 
assess whether the driving of turbulence by supernovae is also important in the outer 
Galactic disk, where the star formation rates are lower, we study the spatial distri- 
bution of molecular cloud (MC) inclinations with respect to the Galactic plane. The 
latter contains important information on the nature of the mechanism of energy in- 
jection into the ISM. We analyze the spatial correlations between the position angles 
(PAs) of a selected sample of MCs (the largest clouds in the catalogue of the outer 
Galaxy published by Heyer et al. 2001). Our results show that when the PAs of the 
clouds are all mapped to values into the [0°, 90°] interval, there is a significant degree 
of spatial correlation between the PAs on spatial scales in the range of 100-800 pc. 
These scales are of the order of the sizes of individual SN shells in low density environ- 
ments such as those prevailing in the outer Galaxy and where the metallicity of the 
ambient gas is of the order of the solar value or smaller. These findings suggest that 
individual SN explosions, occurring in the outer regions of the Galaxy and in likewise 
spiral galaxies, albeit at lower rates, continue to play an important role in shaping the 
structure and dynamics of the ISM in those regions. The SN explosions we postulate 
here are likely associated with the existence of young stellar clusters in the far outer 
regions of the Galaxy and the UV emission and low levels of star formation observed 
with the GALEX satellite in the outer regions of local galaxies. 
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1 INTRODUCTION 

The last few years have witnessed increased efforts in study- 
ing the nature and effects of the different energy injection 
mechanisms into the interstellar medium (ISM) of galax- 
ies. Understanding the nature of the turbulence drivers and 
their associated energy injection rate has important conse- 
quences in improving our understanding of the dissipation 
of turbulence in the ISM and its impact on star formation. 
A few recent studies pointed out that kinetic energy might 
be injected on large scale (kpc scales) into the ISM of dwarf 
irregular galaxies. Stanimirovic & Lazarian (2001) analyzed 
the kinetic energy power spectrum in the Small Magellanic 
Cloud and did not observe any indication of energy injection 
up to the largest considered scale (~ 4) kpc. Dib & Burkert 
(2004,2005) analyzed the HI gas morphology in Holmberg II 
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(Ho II) and a series of numerical simulations of driven tur- 
bulence of the large scale ISM that include, cooling, heating 
and the gas self-gravity. They performed a comparison be- 
tween the observations of Ho II and the simulated HI maps 
by measuring the autocorrelation length on different scales 
and concluded that turbulence is injected into the ISM of 
the dwarf irregular, low star forming galaxy Ho II, on a scale 
of ~ 6 kpc. Similar conclusions on the existence of a large 
scale driving mechanism in the dwarf irregular galaxy DDO 
210 have been reached by Begum et al. (2006). Brunt et al. 
(2009) recently showed that the size-velocity dispersion re- 
lations of nearby Galactic molecular clouds (MCs) can only 
be replicated by models in which turbulence is driven on 
large scales (i.e., larger the sizes of the clouds themselves). 
Koda et al. (2006) measured the orientations of MCs in the 
Galactic Ring region and concluded that the absence of a 
dominant number of clouds oriented perpendicularly with 
respect to the Galactic plane is an indication that turbu- 
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Figure 1. Distribution of the axis ratios for of all molecular 
clouds in the Outer Galaxy molecular cloud catalogue of Heyer 
et al. (1998, 2001) (full line). I max is the size of the cloud along 
its major axis and ,l m i n , the size of the cloud along its minor 
axis. The dashed line corresponds to the sample of clouds with 
Imin ^ 3 pc, l max > 6 pc, and l m ax/l m in > 2. Note that in the 
HCS catalogue l m ax and l m i n fall in the range [0.01,200.7] pc, 
and [0.01,43.5] pc, respectively. 



lence in the Galactic ISM must be injected on scales larger 
than those associated with feedback from massive stars. The 
latter conclusion is based on the fact that supernova (SN) 
explosions occurring in the Galactic disk will generate foun- 
tain like structures and a large fraction of MCs that are per- 
pendicular to the Galactic plane. However, to date, there is 
no complete theoretical/numerical model that includes suffi- 
cient physics and that possesses enough resolution and that 
is able to make predictions about the real fractions of the ori- 
entations of MCs in a Milky Way-like Galaxy. In the absence 
of such information, the conclusions of Koda et al. (2006) 
still require the quantitative argument that would definitely 
confirm or exlude stellar feedback as being responsible for 
the orientations of Galactic MCs. 

The question of what causes the inclinations of MCs 
with respect to the galactic plane in the Galaxy and in 
other galaxies is intimately related to the question of the 
nature of the dominant energy injection mechanism into the 
ISM. Numerical simulations and theoretical studies show 
that feedback from massive stars in the form of stellar winds 
(Krumholz et al. 2006) or SN explosions (Parous et al. 1990, 
Korpi et al. 1999; de Avillez & Breitschwerdt 2005; Slyz et al. 
2005; Joung & Mac Low 2006; Dib et al. 2006; Booth & The- 
uns 2007; Shetty & Ostriker 2008) is able to explain the ob- 
served velocity dispersions of the HI gas in the inner and in- 
termediate regions of galactic disks and a fraction of the ob- 
served velocity dispersion in the outer regions. SN feedback 
also produces substantial numbers of MCs with non-zero 
inclinations with respect to the galactic plane. Using a feed- 
back efficiency of 0.25 (each SN injects E SN = 0.25 x 10 51 
erg) , Dib et al. (2006) showed that feedback from SNe explo- 
sions can maintain a velocity dispersion of 3 — 5 km s _1 of 
the HI gas for SN rates such as those prevailing in the outer 
region of galactic disks; i.e. for SN rates as low as 0.01 times 
the Galactic SN type II rate. Dib et al. (2006) attributed the 
differences between the measured velocity dispersions of the 



HI gas for such low rates and the observed values (~ 5 — 10 
km s~\ Dickey et al. 1990; Dib et al. 2006; Tamburro et al. 
2009) in the outer galactic regions to a number of poten- 
tial effects. The discrepancy can result from the unknown 
temperatures of the HI gas in the outer regions of galax- 
ies and our ignorance about the true thermal component in 
the observed line width, an underestimate of the SN effi- 
ciency, and/or can be due to the presence of other driving 
mechanisms in the those regions, such as the magnetorota- 
tional instability (Selwood & Balbus 1999; Kim et al. 2003; 
Dziourkevitch et al. 2004; Piontek & Ostriker 2007), large 
scale gravitational instabilities, non-axisymmetric perturba- 
tions, and cloud-cloud collisions (Wada et al. 2002; Dib & 
Burkert 2005; Li et al. 2006; Tasker & Tan 2008; Agertz et al. 
2009), tidal interactions and ram pressure effects (Bureau et 
al. 2004); collisions of high velocity clouds with the galactic 
disk (Tenorio-Tagle et al. 1987; Santillan et al. 2007; Baek 
et al. 2008) and galactic spiral shocks (Kim & Ostriker 2006; 
Bonnell et al. 2006; Kim et al. 2008; Dobbs et al. 2008). 

In the present paper, we complement the work of Koda 
et al. (2006) by analyzing the position angles of MCs (i.e., 
their inclinations with respect to the Galactic plane, PA) in 
the outer Galactic disk using the Heyer, Carpenter & Snell 
molecular cloud catalogue (HCS, Heyer et al. 2001). As a 
quantitative test to measure the signature of the turbulence 
driver(s) in the Outer Galaxy, we measure the spatial corre- 
lations between MCs with the same PAs with respect to the 
Galactic plane. In particular, supernova feedback will gen- 
erate a network of MCs of various inclinations. If the PAs 
of the clouds are mirrored to values between 0° and 90°, 
the spatial scales on which they might be correlated can be 
compared to the expected sizes of SN remnants in the outer 
Galactic disk. 

Another possible approach would be to use shell-like 
structures as a direct tracer of feedback from massive stars. 
Shells observed in the HI 21 cm line are not an automatic 
tracer of star formation activity. Although some HI shells are 
generated by feedback from massive stars, Blitz et al. (2007) 
discussed the fact that whereas Giant Molecular Clouds are 
found to be well correlated with high density regions in the 
HI (see also Margulis et al. 1988), the opposite is not true 
and many bright HI filaments are found without molecular 
gas. Dib & Burkert (2005) and other groups have shown that 
HI shell-like structures that are not associated with stellar 
feedback can also be generated by large scale turbulence and 
gas instabilities. In the SMC, Hatzidimitriou et al. (2005) 
found that 59 of the 509 HI shell they have catalogued show 
no signs of being associated to any form of star formation 
or stellar feedback activities. In Holmberg II, 86 percent of 
the HI holes do not show signs of being associated with stel- 
lar feedback (Rhode et al. 1999). This argument does not 
apply to CO shells, as it is very likely that large and com- 
plete CO shells (> 70 pc are associated with feedback from 
massive stars. However, CO shells that would form through 
feedback by massive stars will generally evolve in a turbu- 
lent and highly inhomogeneous medium and may lose their 
spherical symmetry under the action of hydrodynamical in- 
stabilities. This would make their recognition as tracers of 
feedback by massive stars very problematic. More direct ev- 
idence for SN activity would be very difficult to observe in 
the outer Galaxy. Assuming that the SN rate in the outer 
Galactic regions is somewhere between 1/10 to 1/100 of the 
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Figure 2. Distribution of the position angles for the entire sample 
of cores of the HCS catalogue (full line) and the sample of selected 
clouds with l min ^ 3 pc, lmax ^ 6 pc, and l m ax/lmin ^ 2 
(dashed line). 



Galactic value, this would result in an SN type II frequency 
of 1 SN per 500-5000 years since the galactic frequency is 
1/(57 years) (Capellaro et al. 1999). It might also be dif- 
ficult to associate any hot gas with past and present local 
SN explosions as the hot gas will be easily lost to the halo. 
One possibility that remains is to to detect remnant stellar 
clusters associated with the shells. In view of the difficulties 
in the usage of shell-like structure algorithms to trace the 
effects of massive stars feedback, we resort to the method 
described above which is to search for spatial correlations in 
the molecular clouds orientations with respect to the Galac- 
tic disk. 

In § [2] and § [3] we describe the MCs sample and the 
analysis method, respectively. The results are presented in 
§[4] and in §0 we conclude. 



2 MOLECULAR CLOUD SAMPLE 

For our study of the orientations of MCs in the outer Galaxy, 
we use the HCS catalogue (Heyer et al. 1998,2001) which is 
based on the FCRAO 12 CO J=l-0 line survey in the outer 
Galaxy and contains a total of 10156 objects. The survey 
covers 336 deg 2 (40° x 8.4°) in the second quadrant of the 
Galaxy with an angular resolution of 45". An advantage of 
studying MCs in the outer regions of the Galaxy is that their 
distances can be uniquely estimated from the observed ra- 
dial velocity and an assumed rotation curve of the Galaxy 
(in this case a flat rotation curve in the outer Galaxy with 
a circular velocity of v C i rc — 220 km s^ 1 ). Individual MCs 
were identified in this survey as being closed topological sur- 
faces in the I — b — Vlsr space at a given threshold of the 
antenna temperature (1.4 K), where I and b are the Galactic 
longitude and latitude, respectively, and Vlsr is the veloc- 
ity in the local standard of rest. The HCS catalogue lists a 
large number of the MCs properties among which are the 
values of the major and minor axis lengths, l m ax and l m i n , 
respectively, the PA which is denned as being the angle be- 
tween the major axis with respect to the positive Galactic 




Figure 3. A schematic, simplified, view of supernova explosions 
occurring in the Galactic disk and which illustrates the various 
possible position angles of molecular clouds generated by the su- 
pernova explosion shock. A position angle, PA > 90° and its 
mirrored value, PAm, are highlighted. 



latitude axis measured clockwise, the kinematical distance 
d, and the Galactic coordinates. 

The distribution of axis ratios is shown in Fig. [1] Most 
clouds are found to have an axis ratio lmax/lmin > 2 (in 
physical units l max and lmin fall in the range [0.01,200.7] pc, 
and [0.01,43.5] pc, respectively). As pointed out by Heyer et 
al. (2001), clouds with a large axis ratio have more reliable 
determinations of their PAs. Therefore, we chose to analyze 
only clouds which have lmax/lmin ^ 2. A second criterion we 
apply in selecting clouds is that they be large with l m i„ 3 
pc and lmax > 6 pc. This is to avoid using small clouds 
whose inclinations with respect to the Galactic plane might 
significantly differ from the orientation of larger structures 
they would be part of if a lower threshold for cloud selec- 
tion were to be used and/or because their actual inclinations 
might differ from those of the larger structures from which 
they may have fragmented. Note that at similar kinemati- 
cal distances, the axis ratios of smaller clouds will be poorly 
determined compared to those of larger clouds because of 
the smaller numbers of pixels used in the calculation of the 
inertial matrix. Applying these two criteria reduces the to- 
tal number of analyzed clouds to 1485 (R sample). Fig. [2] 
displays the distribution of PAs for all clouds (full line), 
and for the clouds of sample R (dashed line). This figure 
shows that MCs in the outer Galaxy, for both distributions, 
do not have a dominant PA value. They are rather evenly 
distributed among all possible values of the PAs between 0° 
and 180°. 

Shell-like structures can be visually observed in the HCS 
catalogue. Figure 1 in Heyer et al. (2001) which displays the 
distributions of all CO clouds selected in the velocity inter- 
val -110 to -20 km s" 1 between I = 106° and I = 142° shows 
the existence of large shells centered at (l,b) ~ (138°, 1°), 
(134° ,0°), (133.8°, 1°), (123° ,-1°), and (111°,0°). These shells 
seem to extend over 1-2 degrees in longitude. However, it 
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is important to note that the CO shells are incomplete 
which makes the usage of a direct shell-finding algorithm 
very problematic. In the HI, there are also shells in the 
same Galactic region that were catalogued by Ehlerova & 
Palous (2005). The latter authors do not show a large field 
of view of the HI gas distribution, but rather a few se- 
lected HI shells. Among them, the shell that is observed 
at (J,6)=(124° - 131° ,0° - 1°) in figure A.l of Ehlerova & 
Palous seems to coincide in position with two of the CO 
shells. Note that the few HI shells that are shown seem to 
be larger, extending over 4-6 degrees in longitude. 



3 ANALYSIS 

The assumption behind the analysis performed in this work 
is that a supernova or multiple supernova explosions will 
lead to compressions of the surrounding gas that would re- 
sult in the formation of MCs which possess non-zero incli- 
nations with respect to the Galactic plane. Numerical sim- 
ulations of supernova or multiple supernova explosions in a 
stratified galactic disk result in mushroom-like shaped cav- 
ities with dense and cold, molecular material at their edges 
(de Avillez & Breitschwerdt 2005). In order to test the su- 
pernova scenario using the PAs, it is necessary to mirror the 
PAs and restrict them to the range [0°,90°] such that clouds 
that lay on both sides of the remnant with respect to the nor- 
mal direction to the Galactic disk will have their concavity 
pointing in the same direction. This will allow the detection 
of any correlations between clouds on the scale of individual 
supernova remnants or of superbubbles. Hence, we calculate 
the autocorrelation function (ACF) for two cloud samples: 
sample R, and sample Rm which contains the same clouds 
as sample R, but where the PAs of the clouds in the range 
[90°, 180°] are mirrored such as to be restricted to the range 
of [0°, 90°] (e.g., the mirror value of 91° is 89°, see Fig. |3J. 
The ACF, on a physical scale L, is calculated on the discrete 
sample of clouds using the usual definition: 



S 2 (L) = 



EJ£Mf) - PA)(PAj(r + L) - PA) 
£ t (JM<(r) - PA)» 



(1) 



where the index i runs over all MCs in the used sample, PA 
is the average PA in the sample and L is a physical lag. 
The ACF is evaluated over a series of lags incremented by 
100 pc. Houlahan & Scalo (1990) pointed out that, when ap- 
plied to extended structures, the ACF can be separated into 
self-terms which are sensitive to their substructure and to 
to cross-terms which are more sensitive to their separation. 
Houlahan & Scalo (1990) showed that the cross-terms can 
be overshadowed in a hierarchical-like structure by the self- 
terms and by other effects such as boundary effects. They 
point out nevertheless that the usage of the ACF should not 
be problematic in the case of point like structures, and that 
in this case the ACF is essentially sensitive to the point- 
point separation. In our case, since we do not apply the 
ACF on every single pixel of the clouds (i.e., we do not 
use their extended gas distribution), we are in the regime 
of point sources since every cloud is only represented by a 
unique information which is the PA. As we check from the 
position of each cloud what are the other clouds that are at 
distances between r and r + L from this cloud, it is necessary 
to know the cloud-cloud distances. The latter are calculated 
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Figure 4. Kincmatical distances to the clouds of the entire 
sample of clouds (full line) in the HCS catalogue (Heyer et al. 
1998, 2001) and the sample of selected clouds with l m in 3 pc, 
' ni a x pc, and l-max /lmin ^ 2 (dashed line). 
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Figure 5. Spatial autocorrelation function of the position angles 
of the selected clouds in the HCS catalogue (i.e., clouds with 

lmin ^? 3 pc, Ijnax 6 pc), and Imax I lmin z£ 2. 



using the Galactic coordinates and the kinematical distance 
to each cloud and the cosine law in spherical trigonometry. 
Fig.|4]displays the kinematical distances to the clouds in the 
HCS catalogue as derived by Heyer et al. (2001) (full line) 
and for those of the R (and Rm) samples (dashed line). The 
cloud-cloud distance, dij, between cloud i and cloud j which 
are located at distances D; and Dj and which have Galactic 
coordinates (h,bi) and (lj,bj), respectively, is given by: 



dij = [D~ + D 



2 DiDjCos(bi)cos(bj)cos(li — lj) 
—2 DiDjsin(bi)sin{bj)] L ^ 2 . 



(2) 



4 RESULTS 

Fig. [5] displays the spatial ACF of the PAs for the clouds 
of sample R. The figure clearly shows that there is no sub- 
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Figure 6. Spatial autocorrelation function (ACF) of the position 
angles mirrored into the range [0,90]° for the selected clouds in 
the HCS catalogue (clouds with l min ^ 3 pc, Imax ^ 6 pc, and 
Imax /lmin 2. The dashed line corresponds to the ACF for 
the same population of clouds but with random position angles 
(built from the average of 10 individual ACFs with 10 random 
seed numbers). 

stantial correlation between the orientations of the PAs for 
the clouds of sample R on any physical scale. On the other 
hand, Fig. [S] displays the spatial correlations of the P^4s for 
the Rm sample. In contrast to the un-mirrored case, the 
mirrored PAs for the selected clouds show non-negligible 
correlations on spatial scales of the order of ~ 100 — 800 pc. 
We have also checked for the consistency of the results by 
taking different permutations of l m in and Imax and of their 
imposed ratio; i.e., values of (l m i„ (j>c),l m ax (pc)) = (2,5), 
(3,5), (3,7) with and without imposing that Imax/lmin ^ 2, 
and found no significant variations in the results for Fig. [S] 
and Fig. |§] These spatial scales of a few hundred parsecs 
are well matched by the sizes of supernova remnants in low 
density environments. Indeed, using analytical calculations, 
McCray & Kafatos (1987) estimated the cooling radius of a 
supernova shell/supershell to be given by: 

i? c = 50C- a9 (iV,i5 5 1 4 no - 6 )p C , (3) 

where C, is the local value of the metallicity, N* is the num- 
ber of stars in the stellar cluster with masses greater than 
7 Mq, E$\ is the energy input of one supernova in units of 
10 51 erg, and no, in cm~ 3 , is the ambient gas density. The 
cooling radius marks the end of the adiabatic phase of the 
expanding shell/supershell and the time at which cooling 
becomes important in the hot interior of the remnant. Note 
that the shell/supershell continues to expand afterwards ac- 
cording to the zero-pressure snowplow law with its radius 
slowly increasing following R(t) = i? c (f/f c ) 1/ ' 4 , where t c is 
the time at which R c has been reached. For an ambient gas 
density no = 0.1 cm -3 , AT* = 1, and a solar metallicity, the 
value of R c is of the order of ~ 200 pc, for N* = 1, ( = 1, 
and no = 0.03 cm -3 , R c ~ 410 pc, whereas for no = 0.03 
cm" 3 , N* = 1, and a lower metallicity value of £ — 0.5, 
the cooling radius is R c ~ 765 pc. In order to compare the 
correlations to the noise, Fig. [5] also displays (dashed line) 
the ACF for the same sample of clouds but with their PAs 
drawn from random values (averaged over 10 iterations). As 



expected in this case, the ACF tends to zero values on all 
spatial scales. 

The similarity between the latter physical scales and 
the physical correlation scales observed in Fig. [5] suggests 
that supernova driving might be playing a considerable role 
in shaping the ISM structure and dynamics even in regions 
of low star formation activity. This is consistent with the 
results of Dib et al. (2006) who showed that supernova feed- 
back is capable of maintaining turbulent velocity dispersions 
for the HI gas of the order of 3 — 5 km s _1 even in regions 
where the star formation rate per unit area is as low as 
10~ 4 Mq yr _1 kpc -2 , corresponding to one hundredth of 
the Galactic SN rate. The existence of spatial correlations 
for the sample of mirrored PAs which we associate to SN 
shells and the absence of such correlations for sample R in- 
dicates that there are no correlations between the sites of 
star formation and that gravitational instability occurs in 
the Galaxy wherever local conditions are favorable. 

Note that, in calculating the ACFs, we have made use 
of the kinematical distances to each of the clouds published 
by Heyer et al. (2001). The kinematic distance in the outer 
Galaxy is an overestimate (e.g., Brand & Blitz 1993) of the 
true distance measured by trigonometric parallax measure- 
ments (e.g., Hachisuka et al. 2006,2009). For example, more 
accurate distance estimates using parallax measurements for 
the W3 region yield a distance of 2.1 kpc (Hachisuka et al. 
2006) which is smaller than the value of ~ 3.5 kpc we have 
used for clouds in that region. The overestimate of the dis- 
tance by the kinematic distance in the outer Galaxy will 
depend on the position of the cloud and on the existence of 
random velocity components and/or velocity offsets at this 
location with respect to a given Galactic rotation curve. We 
have calculated simple models to quantify the kinematic dis- 
tance uncertainties when using the flat rotation curve (with 
a circular velocity of v C i rc — 220 km s _1 similar to the value 
adopted in Heyer et al. 2001) for objects within the longi- 
tude range studied in this paper. The method proceeds as 
follow: We first calculate a distance assuming pure circular 
motions for a range of expected velocities and longitudes 
in the ranges [-5,-100] km s _1 and [100°, 140°], respectively. 
In a second step, we add a velocity offset (Aw) and/or a 
random component with a dispersion of u v and recompute 
the distance. The procedure is repeated 8192 times in order 
to sample the random component of each longitude and ve- 
locity. With that, we can derive the root mean square, <jd, 
of the newly calculated distance. Fig[7] shows the fractional 
distance error od/D for three sets of a v and A„ as a func- 
tion of the velocity for various values of the longitude I. If 
there are no streaming motions (A„ = 0) but only a ran- 
dom velocity component of dispersion a v = 5 km s _1 (left 
plot in Fig[7J), the differences between the kinematic distance 
and the true distance are of the order of < 10 percent for 
most velocities and get larger for nearby clouds. Such errors 
could be propagated into the calculation of the inter-cloud 
distances using Eq[2] on a two-by-two basis. However, the 
overall uncertainty on the ACF is not expected to be very 
large and might stretch, for a given physical scale, the x-axis 
in Fig[5]and Fig[6]by a factor of 10 percent in one direction 
or the other. Streaming motions introduce larger but sys- 
tematic uncertainties which in fact are overestimates of the 
true distance (models with Av — 15 km s _1 , Fig. [7] middle 
and right). The systematic overestimates can be larger than 
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the true distance by factors of up to a few. Since most of the 
clouds in the sample have velocities v < —25 km s -1 , ctd/D 
is expected to be < 1. To correct for these distance overesti- 
mates, the x-axis in Fig. [5] and Fig. [6] should be contracted 
by a factor of \{<td / D) + 1] ~ 2. This would cause the mir- 
rored PAs to be correlated on scales between ~ 50 — 400 pc. 
Nevertheless, this would not affect our conclusions. 

In conclusion, our results suggest that SNe do occur in 
the outer regions of the Milky Way (MW). That stars form 
at large galacto-centric radii in MW like galaxies is sup- 
ported by the discovery of young stellar clusters at larger 
Galactocentric distances (e.g., Carramaza et al. 2008; Yun 
et al. 2009) and UV observations by the GALEX satellite. 
The latter show evidence for extended UV emission presum- 
ably related to ongoing star formation and feedback from 
massive stars, at low rates, in those regions (Gil de Paz et 
al. 2007). This star formation extends beyond the region of 
the canonical star formation threshold (Martin & Kennicutt 
2001). Observationally, H a emission, which is indicative of 
star formation activity is not always associated to the UV 
emission (Gil de Paz et al. 2007). In fact, Boissier et al. 
(2007) argue that because of the extremely low star forma- 
tion rates, the most massive part of the initial mass function 
is only sparsely sampled and not all star forming complexes 
host massive stars. SN explosions are at least expected to 
be associated to a fraction of the star forming sites in the 
outer galactic regions. Using the Leiden-Dwingeloo HI sur- 
vey, Ehlerova & Palous (2005) analyzed the properties and 
radial distribution of HI shells in the second quadrant of 
the outer Milky Way. They found that the radial decrease 
in the surface density of shells in the outer Galaxy is ex- 
ponential, with a radial length scale of ~ 3 kpc and which 
is comparable to the length scale of the stellar distribution. 
The similarity between the length scales of stars and HI 
shells suggests that at least some shells must be connected 
to stars and to the effects of stellar feedback. Elmegreen & 
Hunter (2006) showed that the radial star formation rates 
in the disks of spiral galaxies decline gradually with radius 
rather than being subjected to a sharp cutoff. This is due 
to the persistence of various gas instabilities in the outer 
regions such as stellar feedback, gravitational instabilities, 
and spiral shocks. 



5 CONCLUSIONS 

We analyze the spatial correlations of the positions angles 
(PAs) of molecular clouds, i.e., their inclination with respect 
to the Galactic plane, for a sample of CO clouds located 
in the outer Galaxy (the CO cloud catalogue of Heyer et 
al. 1998,2001). When the PAs of the clouds that fall in the 
range [90°, 180°] are mirrored to values in the range [0°,90°], 
there is a significant degree of spatial correlations between 
the PAs of the largest molecular clouds on spatial scales 
in the range of ~ 100-800 pc. These physical scales are of 
the same order as the expected physical sizes of individual 
supernova (SN) shells in low density environments such as 
those prevailing in the outer Galaxy and where the metal- 
licity of the ambient gas is of the order of the solar value or 
smaller. Our results suggest that individual SN explosions 
continue to occur in the outer regions of the Galaxy, even at 
lower rates, and may play an important role in shaping the 



structure and dynamics of the interstellar medium in those 
regions. 
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Figure 7. Fractional distance errors for clouds in the outer Galaxy (second quadrant) as a function of their velocity and longitude 
assuming that (left) clouds have a random velocity component a v = 5 km s -1 and no streaming motions Av = around the fiat circular 
velocity v c i rc = 220 km s — 1 , (middle) clouds have streaming motions Ad = 15 km s — 1 and no random motions (<r„ = 0) around the 
flat circular velocity v c i rc = 220 km s _1 , and (right) clouds have a random velocity component a v = 5 km s -1 and streaming motions 
with Av = 15 km s — 1 around the flat circular velocity v c i rc = 220 km s — 1 . The solid, dotted, dashed, dot-dash and triple dot-dash lines 
correspond to the longitude values of I = 100°, 110°, 120°, 130°, 140°, respectively. Note that in the existence of streaming gas motions 
(middle and right figures), the kinematic distance corresponds always to an overestimate. 
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